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Abstract 



The cubic lattice symmetry of ferromagnetic homogeneously disordered al- 
loys is broken when a compositional modulation is imposed. This can have a 
profound influence on the magnetocrystalline anisotropy energy (MAE). We 
describe our ah initio theory of this effect and use the framework of concentra- 
tion waves with the electronic structure described within the spin-polarised 
relativistic Korringa-Kohn-Rostoker coherent-potential approximation. We 
find that ordering produces a 2 order of magnitude increase in the MAE as 
well as altering the equilibrium direction of magnetisation. Using the same 
theoretical framework we also examine directional compositional order pro- 
duced by magnetic annealing with an explicit study of permalloy. 



PACS numbers: 75.30.Gw, 75.50.Cc, 75.60.Nt, 75.50.Ss 
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Magnetocrystalline anisotropy (MCA) of ferromagnetic materials containing transition 
metals has become the subject of intensive theoretical and experimental research because 
of the technological implications for high-density magneto-optical storage media P]||. Po- 
tential materials for these applications need to exhibit a substantial perpendicular magnetic 
anisotropy (PMA), which is evidently due to an intrinsic magnetic anisotropy in the crys- 
tal lattice strong enough to overcome the extrinsic macroscopic shape anisotropy favoring 
an in-plane magnetisation. Whereas in ultrathin films and multilayers the PMA is due to 
surface and interface effects respectively, in thick films of transition metal alloys 
it is an intrinsically bulk magnetic property which leads to PMA. Such systems are par- 
ticularly interesting for magneto-optic recording because in addition to PMA they exhibit 
large Kerr effect signals as well as being chemically stable and easy to manufacture. There 
has been much effort directed towards an understanding of the mechanism of MCA from a 
first-principles electronic structure point of view to aid future magnetic material design but 
since MCA arises from spin-orbit coupling, which is essentially a relativistic effect || this 
means that a fully relativistic electronic structure framework is desirable. 

Empirically compositional structure is found to have a profound influence on both the 
magnitude of the magnetocrystalline anisotropy energy (MAE) as well as the equilibrium 
or easy magnetisation direction. Compositional order lowers the lattice symmetry of the 
homogeneously disordered alloy and this enhances its MAE. For example, the measured 
MAE (~ 130/ieV) of ordered CoPt alloy (CuAu or Ll type) [|7| is some two orders of 
magnitude larger than that of its disordered counterpart (~3 /xeV), and with a different easy 
axis. In this letter we provide the first quantitative calculations of the effect of compositional 
order on MAE via a study of CoPt. 

We also note that thick fcc-Co c Pti_ c @-[K| and fcc-Co c Pd!_ c [g films exhibit PMA of a 
size comparable with that of multilayers. Such a large PMA is quite unexpected for systems 
which should have effectively bulk cubic structure. Now whereas tensile strain and large 
negative magnetostriction coefficients may account for some of this in the Co c Pdi_ c films, 
it cannot for the thick Co c Pti_ c ones Rather recent experiments have confirmed that 



although Co c Pti_ c films seem to be nearly homogeneously disordered there are actually more 
Co-Co nearest neighbors in-plane and very few out-of-plane [y,|l^,|n|, i.e. there is presence of 
some atomic short-range order (ASRO). The existence of this compositional order produces 
internal interfaces analogous to Co/Pt multilayers and this has been suggested as a likely 
cause of the strong PMA in fcc-Co c Pti_ c films [|K|[12|- This suggestion is strengthened by 
the fact that films grown at higher temperatures, in which enhanced bulk diffusion tends 



to destroy the in-plane compositional order, show no PMA. Very recently Kamp et al [14 



have shown from magnetic circular x-ray dichroism measurements that chemical ordering 
is also responsible for the enhanced MAE in thick Feo.sPdo.s films. These observations 
clearly underline a correlation between MCA and ASRO. In this letter we provide, for 
the first time, a detailed theoretical analysis based on first-principles calculations which 
shows how chemical order significantly influences MCA and compare the MAE of ordered 
and disordered CoPt alloys. Moreover, we use the same approach to model other hitherto 
unfabricated structures. This is pertinent now that it is possible to tailor compositionally 
modulated films to obtain better magneto-optic recording characteristics [|T5[|. 

Upon ordering CoPt undergoes a modest tetragonal lattice distortion (c/a = 0.98) which 
also lowers the symmetry. From magnetostriction data however we can estimate this to 
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contribute only 15% of the MAE. The measured magnetostriction coefficient, A, of pure 
fcc-Co (5 x 10~ 5 ) is of the same size as that of CoPt (-4 x 10~ 5 ) although with different sign. 
From 'first-principles' electronic structure calculations Wu et al |16|] have shown that the 



rate of change of the MAE with lattice strain is proportional to A and thus in fcc-Co and 



CoPt should be roughly of the same magnitude. Thus from Wu et a/'s [[16|] calculation of 
the MAE of Co for a range of c/a ratios we estimate that a 2% tetragonalization will change 
the MAE only by about 20/ieV. Consequently we conclude that it is the compositional order 
that is primarily responsible for the large MAE of CoPt. 

In previous work, we presented a theory of MCA of disordered bulk cubic alloys 
fT7| , [T8|| within the framework of spin-polarised relativistic Korringa-Kohn-Rostoker coherent- 
potential approximation (SPR-KKR-CPA) In this letter we set up a new framework 



to investigate the effects of compositional order, both short and long-ranged. Any com- 
positionally modulated alloy can be specified by site-dependent concentrations {q} which 
themselves can be written as a superposition of static concentration waves (CWs) PD| , i.e., 



1 q 

with wave-vectors q and amplitudes c q . Usually only a few CWs are needed to describe a 
particular ordered structure. For example, the CuAu-like Ll structure (Fig. 1) is set up by 
a single CW with c q = | and q = (001), and the [lll]-layered CuPt-like L\\ structure is set 
up by a CW with c q = | and q = {\\\) (q is in units of — , a being the lattice parameter). 

The grand-potential for the interacting electrons in an alloy with composition {cj} and 
magnetised along the direction e x at a finite temperature T is given by pi, 22 . 



/oo 
def(e, u)N({d},e; e x ) 
-oo 

+ Q DC ({ci}; ei), 

where, v is the chemical potential, Z the total valence charge, f{s,i>) the Fermi factor, 
N({ci}, e; ei) the integrated electronic density of states, and QDc({ci}',ei) the 'double- 
counting' correction. The MAE of the inhomogeneous alloy can be characterised by the 
difference 

K({c l }) = Q({c t };e 1 )-Q({c l };e 2 ). 

Assuming that flDc{{ c i}', e ) is generally unaffected by the change in the magnetisation 
direction, we get, 

/oo 
def(e, Vl ) [N({*}, e; e x ) - N({ Cl }, e; e 2 )] 
-oo 

+ 0{y x - v 2 f. 

Note that the correction due to the change in the chemical potential (from v\ to V2) with 
the magnetisation direction is of second order in iv\ — z/2), and can be shown to be very 
small compared to the first term ||17|| . We now expand K({ci}) around Kcpa{c), the MAE 



of the homogeneously disordered alloy A C B 



l-c, 
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K({ Ci }) = K CPA {c) +£ 



dK({ Ci }) 



5cj 



1 d 2 K({ Ci }) 
^ dcjdck 



{Ci=C> 

5c,-5c A + 0(<5c) 3 . 



(1) 



{ci=c} 



The second term in Eq. (P vanishes if the number of A and B atoms in the alloy is to 
be conserved. Now taking the Fourier transform, we get the MAE of the compositionally 
modulated alloy with wave-vector q, 



K{q) = K CPA (c) + -\ C(l \ 



^(qjeO-SOOCqje,,) 



(2) 



and so, for example, the MAE of the CuAu-type Lin ordered alloy is obtained by choosing 
q = (001) and c q = ~. Here 5^(q; e) are the Fourier transforms of the so-called direct cor- 
relation functions |HJ] which determine the ASRO parameter a(q) in the disordered phase, 
a(q) = /3c(l — c)/[l — (3c(l — c)S^ 2 '(q)], {(3 = l/kgT, fee being the Boltzmann constant). 
These have been calculated for many alloys, both non-magnetic and ferromagnetic, in which 
up to now relativistic effects were largely ignored and were compared with diffuse X-ray and 
neutron scattering data [pl],Q. S^(q) is given by, 



Im 
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/oo 
def(e,is) Yl (X A -X 
OO T r T T 



L1L2L3L4, 



(q)A i3L4 (q) 



where, 



^L 1 L 2 (<1) — (X A - X B ) LlL2 



l^uwmu (q)^f 6 L 2 ^l 3 l 4 (q) 



q) = t^- f dkT L5L3 (k + q)T L4L6 (k) 
vbz J 



no 



no 



T L 5 L 3 T L 4 L 6 ^ 



(3) 



and X A( - BS) = [(tJu B ) — tc 1 ) 1 + t 00 ] 1 - Here, r 00 is the site-diagonal path-operator matrix, 
£a(b) and t c are the t-matrices for electronic scattering from sites occupied by A(B) atoms 
and the CPA effective potentials respectively, r(k) = [t^ 1 — ^(k)] -1 , and g(k) is the KKR 
structure constants matrix The spinodal ordering temperature below which the alloy 
orders into a structure characterised by concentration wave- vector q ma x is given by pljp2| . 
T c = c(l — c)S^ (q^max', e)//uB, where q max is the value at which 5"^(q;e) is maximal. 

As with all MAE calculations, given the sizes of the energies involved, numerical com- 
putation of K(q) needs to be done very carefully. The energy integration is done using a 
complex contour described elsewhere [17]. The Brillouin zone (BZ) integration is done using 
the adaptive grid method [24"] which has been found to be very efficient and accurate. In 



this method one can preset the level of accuracy of the integration by supplying an error 
parameter e. The integration in Eq. @ is done with e = 10 -6 which means that S^^q; e) 
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which are of the order of 0.1 eV are accurate up to 0.1 fieV. To achieve such level of ac- 
curacy, we had to sample a large number of k-points in the BZ. Also, owing to the form 
of the integrand in Eq. (|3[) the integration has to be done using the full BZ. Typically, in 
our calculations, we needed around 35,000 k-points for points on the energy contour 0.5 Ry 
above the real axis. When the energy was 0.0001 Ry above the real axis (and that is the 
closest point to the real axis, we need) we required as many as 3 million k-points for the same 
level of accuracy. Furthermore, we have calculated S^(q; ei) and S^(q; simultaneously 
ensuring that they are calculated on the same grid, and thus cancelling out the systematic 
errors if any. Therefore, we claim that the values of K(q) are accurate to within 0.1 /ieV. 

We summarise the results for fcc-Co .5Pt .5 alloy in Tables |and [TT]. We note that, S'^(q) 
is maximum for the LIq structure, implying that the alloy would order into this structure at 
1360 K as it is cooled from high temperature in good agreement with experiment (ordering 
temperature of 1000 K 0). From Table [Ij we note that for q = (001) and q = (100) which 
correspond to CuAu-like Ll ordered structure, with Co and Pt layers stacked along the 
[001] and [100] directions respectively, the direction of spontaneous magnetisation is along 
the [001] and [100] directions respectively in excellent agreement with experiment. Also, the 
MAE (58.6 /ieV) is comparable to the experimental value (~ 130/ieV 0). 

In order to probe the relationship between the compositional structure and MAE further, 
we also performed calculations for some hypothetical structures. These are also summarised 
in Table |TJ. The case of q = (§§|) corresponding to the CuPt-type Ll\ ordered struc- 
ture with Co and Pt layers perpendicular to the [111] direction of the crystal produces a 
spontaneous magnetisation along the [111] direction of the crystal. This may be close to 
the structure of the thick [111] oriented disordered CoPt alloy films exhibiting PMA [p[-|T0| 
which is attributed to the existence of internal interfaces [1^|, analogous to Co and Pt lay- 
ers along the [111] direction. Our result is clearly consistent with these observations. The 
point is that the magnetic anisotropy (152/zeV) intrinsic to this structure is nearly 3 times 
larger than that of Ll structure. Indeed, we predict that a [lll]-oriented film will exhibit 
a markedly stronger PMA than that of a [100]-oriented film. 

The structure set up by CWs with q = (10|) and (0l|) is also a [001] -oriented layered 
structure, but the layers are not alternately pure Co and Pt planes, rather they are layers 
of ordered Co and Pt (Fig. 1). Even in this case, we note that the magnetisation is 
perpendicular to the layered structure, and the magnitude of MAE is large. Similarly, for 
q = (|01) and (|l0) where the planes are stacked along the [100] direction the magnetisation 
is also along the [100] direction. 

Evidently, the spontaneous magnetisation always tends to align itself perpendicular to 
any layering in the structure. In addition, the magnitude of MAE depends strongly on the 
symmetry of the system, i.e. it increases when the symmetry is lowered. The cubic symmetry 
of the homogeneously disordered alloy quenches the orbital magnetic moment. In the LIq 
ordered alloy the symmetry is lower and in the Lli layered structure it is lower still, thus 
increasing the MAE in each case. On detailed examination of the electronic structure of the 
disordered alloy around the Fermi energy we find that the large values of K(q) near the BZ 
boundary arise from van Hove singularities of the Bloch spectral function [p3fl . The number, 
location, and occupation of these depend on the magnetisation direction and produce a large 
contribution to the difference in the convolution integrals (Eqs. (Q) and (|3|)). 

Our theory can also be used to produce the first quantitative description of the phe- 
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nomenon of magnetic annealing. Here a soft magnetic material develops directional chemi- 



cal order when annealed in a magnetic field ||25|| . We demonstrate this effect for Nio.75Feo.25 
(permalloy) and Table [TTI| is a summary of the results. We calculate (q, e) for permal- 
loy in an applied magnetic field of strength 600 Oe (same as used in the experiment p5| ) 
orientated along e = [001], [111], and [100] directions of the crystal. We find that when the 
magnetic field is along [001] (column 2) S^(q) is maximum for q=(001) confirming that 
ordering is favored along the direction of applied field. Similar is the case when the magnetic 
field is along [100]. However, when the applied field is along [111] (column 6) all the three 
orderings, namely, (100), (010), and (001) are favored equally. Thus, in this case, we will 
get a CuAu 3 -type Ll 2 ordering. The calculated transition temperature 721 K is in good 
agreement with the experimental value of 820 K p6 |. Noting that the measured intensity in 
a scattering experiment is proportional to the ASRO parameter a(q), we estimate that for 
an alloy cooled in a magnetic field along the [001] direction the superlattice spot at q=(001) 
will be 20% more intense than that at q=(100) at a temperature 1 K above the transition 
temperature. 

In conclusion, we have presented a fully relativistic electronic structure scheme to study 
the MCA of alloys and its dependence upon compositional structure. We applied this theory 
to fcc-Coo.sPto.s and found that when the system is cooled it tends to order into Ll layered- 
ordered structure and that the spontaneous magnetisation tends to align itself perpendicular 
to the layer stacking. These observations are in complete accord with experiment. We also 
found that if the layers are stacked along the [111] direction then the MCA becomes larger 
which may be the case in [lll]-textured films. Within the same framework we have also 
been able to explain the appearance of directional order in Nio.75Feo.25 when it is annealed 
in a magnetic field. 

We thank B.L. Gyorffy for many helpful discussions. This research is supported by the 
Engineering and Physical Sciences Research Council (UK), National Science Foundation 
(USA), and the TMR Network on "Electronic structure calculation of materials properties 
and processes for industry and basic sciences" . 
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TABLES 



TABLE I. Direct correlation function S'^^q; [001]) for different q-vectors for Coo.sPto.s alloy. 







s (2) (q; [ooi]) 


T c 


q 


Structure 


(eV) 


(K) 


(000) 


Clustering 


-1.51 




(001) 


Ll 


0.47 


1360 


(100) 


Ll 


0.47 


1360 



(iii) Lli 0.29 

(10*) 0.19 



TABLE II. Magnetocrystalline anisotropy energy K(q) for several compositionally modulated 



CoPt alloys characterised by different q-vectors (the respective ordered structures are shown in 
Fig. 1). The reference system has the magnetisation along the [001] direction (i.e. ei = [001]). 



q 


e 2 = [111] 
K(d) 
(/xeV) 


e 2 = [100] 
K(d) 
(/zeV) 


Easy Axis 


(001) 
(100) 

fill) 

V2 2 2> 

(10*) 
(±01) 


-58.6 
39.6 
152.0 
-158.7 
85.0 


-105.6 
105.9 
0.0 

-236.5 
236.3 


[001] 
[100] 

[111] 

[001] 
[100] 



TABLE III. Direct correlation function S^ 2 \q;e) for different q-vectors for Nio.75Feo.25 alloy. 



q 


e=[001] 

s (2) (q) 


T c 


e=[100] 
S (2) (q) 


T c 


e=[lll] 

S (2) (q) 


T c 




(meV) 


K 


(meV) 


K 


(meV) 


(K) 


(100) 


330.880 




331.158 


720.9 


330.973 


720.5 


(010) 


330.880 




330.880 




330.973 


720.5 


(001) 


331.158 


720.9 


330.880 




330.973 


720.5 


(111) 

V2 2 2> 


-99.291 




-99.291 




-99.475 




(10*) 


99.284 




99.096 




99.222 




(2-01) 


99.096 




99.285 




99.222 
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FIGURES 

FIG. 1. Some ordered structures and their representative wave-vectors. q=(001) generates the 
CuAu-type structure with ordering along the [001] direction. f\={\\\) generates the CuPt-type 
structure. q=(10^) generates a layered structure with planes of an ordered structure of A and B 
atoms stacked along the [001] direction. 
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